LU TP 01-22 
June 2001 

The Lund Fragmentation Process for a 
Multi-gluon String According to the Area Law 

Bo AnderssonJ] 

| Sandipan Mohantyn 

■ Fredrik SoderbergjJ 

Cn| I Department of Theoretical Physics, Lund University, 

£ ■ Solvegatan 14A, S-223 62 Lund, Sweden 

Abstract 

\Q ' The Lund Area Law describes the probability for the production of a set of 

colourless hadrons from an initial set of partons, in the Lund string fragmen- 
tation model. It was derived from classical probability concepts but has later 
been interpreted as the result of gauge invariance in terms of the Wilson gauge 
loop integrals. In this paper we will present a general method to implement the 
qq , Area Law for a multi-gluon string state. In this case the world surface of the 

massless relativistic string is a geometrically bent (1 + l)-dimensional surface 
embedded in the (1 + 3)-dimensional Minkowski space. The partonic states are 
in general given by a perturbative QCD cascade and are consequently defined 
only down to a cutoff in the energy momentum fluctuations. We will show 
that our method defines the states down to the hadronic mass scale inside an 
analytically calculable scenario. 

We will then show that there is a differential version of our process which 
Oh, is closely related to the generalised rapidity range X, which has been used as a 

measure on the partonic states. We identify A as the area spanned between the 
directrix curve (the curve given by the parton energy momentum vectors laid 
out in colour order, which determines the string surface) and the average curve 
(to be called the P-curve) of the stochastic A-curves (curves obtained when the 
5^ ' hadronic energy-momentum vectors are laid out in rank order) . Finally we show 

that from the X -curve corresponding to a particular stochastic fragmentation 
situation it is possible to reproduce the directrix curve (up to one starting vector 
and a set of sign choices, one for each hadron). This relationship provides an 
analytical formulation of the notion of Parton-Hadron Duality. The whole effort 
is made in order to get a new handle to treat the transition region between where 
we expect perturbative QCD to work and where the hadronic features become 
noticeable. 
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Figure 1: A high-energy string breakup. 

1 Introduction 

The Lund String Fragmentation Model was developed many years ago, |ij ,Q, 
and as implemented in the well-known Monte Carlo simulation program JET- 
SET, 0], it has been very successful in reproducing experimental data from high 
energy multi-particle processes. 

The model is based upon a few general assumptions: (i) the final state 
particles stem from the breakup of a string-like force field spanned between the 
coloured constituents, (ii) there is causality and Lorentz invariance and (iii) the 
production of the particles can be described in terms of a stochastic process 
which obeys a saturation assumption. We have, in a recent paper ||, re-derived 
the major result for the (l + l)-dimensional model, which is applicable for events 
with a quark (q, a colour-3) and an antiquark (q, a colour-3) at the endpoints 
of the string but with no interior gluonic (g, colour-8) excitations. The result is 
that the (non-normalised) probability for the production of an n-particle final 
state of hadrons with energy momenta {pj} and masses {rrij} is given by the 
Lund Area Law: 

n n 

dP n ({Pjh Ptot) = J] N,d 2 V ^ V ) - mtyC^Pj - Ptot) exp(-bA) (1) 
i=i i=i 

where A is the area spanned by the string "before" the breakup, cf. Fig.[j], 
Ptot is the total energy- momentum of the state and {Nj} and b are parameters 
related to the density of hadronic states and the breakup properties of the string 
field, respectively. 

The result in Eq. (Q) is evidently similar to a quantum mechanical transition 
probability. It is the final state phase space multiplied by a squared matrix 
element, which in this case would correspond to the negative area exponential. 
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In la] we have shown that it is possible to "diagonalise" the model, i.e. to 
express the result in Eq. (|l|) as a product of (diagonal) transition operators 
(in quantum mechanical language it would correspond to density operators). It 
turns out that in this way the dynamics can be described in terms of (1 + 1)- 
dimensional (space-like) harmonic oscillators. 

Shortly after the original derivation of the Area Law, [El, Sjostrand, Q, 
provided an implementation of the model applicable also for multi-gluon states, 
i.e. when the string surface is no longer flat but geometrically bent due to 
the internal excitations. Sjostrand's method is to project the positions of the 
breakup points (the vertices) from the (flat) (1 + l)-dimensional model as given 
by Eq. (|l]) onto the surface of the bent string. The projection is done so that the 
proper times of the vertices and the squared masses of the particles produced 
between them are the same. Unfortunately this method does not fulfil the Area 
Law on the bent surface because it is a geometrical fact that the areas "below 
the vertices" are not invariant under such a projection from a flat to a bent 
surface. 

Although the Area Law is not fulfilled on an event to event basis by the 
method in j(| we will show that it is fulfilled in an average sense, i.e. the 
predicted inclusive distributions are little affected by the differences. It is well- 
known that the experimental results for these distributions are well described 
by JETSET even up to the largest available energies of today. 

The intention of this note is to implement another method for particle pro- 
duction in multi-gluon states which fulfils the Area Law at every single step in 
the production process. We will find that it is necessary to tackle a set of prob- 
lems in the definition of the states which we apply the process to. We note that 
the states defined by perturbation theory are resolved only to the scale of some 
virtuality cutoff. We will find that our method provides a set of excitations on 
the hadronic mass scale, in the string field. We will investigate the properties 
of these "soft hadronisation gluons" in future work. 

The states of the massless relativistic string fulfil a minimum principle, i.e. 
the surface spanned in spacetime by the string during its motion is a minimal 
surface. This means on the one hand, that the states should be stable against 
small-scale variations and on the other hand that the surface is fully determined 
by the boundary curve. In this case the boundary curve corresponds to the orbit 
of one of the endpoints, conventionally chosen as the g-endpoint. Therefore, the 
process we are going to define is a process along this curve, to be called the 
directrix curve, which is completely defined by the perturbative cascade. In this 
paper we will treat the partons as massless, although both the process and the 
directrix curve can be defined for a general case with massive quarks. 

One property which can be derived from Eq. (|l|) is that the average decay 
region is a typical hyperbola. On the average, the final state hadrons in our pro- 
cess will be produced in the same way, albeit this time along a set of connected 
hyperbolae. In Q we have defined such an average curve and we will, in this 
paper, call it the A'-curve. Just as a simple hyperbola has a length proportional 
to the hyperbolic angle that it spans (this corresponds to the available rapidity 
range along the mean decay region, in a two-jet system of hadrons) the A"-curve 
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has a length corresponding to a generalised rapidity range, usually called the A 
measure, ||, ||]. The A"-curve is defined in terms of differential equations and 
we will show the close relationship between the A"-curve and our process in the 
limit of a vanishing hadron mass. 

There are several reasons to undertake this investigation. One is to compare 
the precise implementation of the Area Law to the approximate process in Q . 
We will do this both in this paper and in future publications. 

Another reason is to get a handle on the general structure of fragmentation, 
in particular to be able to treat also the multi-gluon fragmentation states by 
the analytical methods introduced in 0]. This is of particular interest for the 
transition region, i.e. the region in between where we expect perturbation theory 
to work and where we know that the non-perturbative fragmentation sets in. 

A final reason is to investigate the stability of the states in QCD under 
fragmentation, i.e. given a multi-gluon state defined according to the rules of 
perturbation theory (with cut-offs as mentioned above) to find out to what 
extent it can be modified so that the observable results after fragmentation 
are still in agreement with the experiments. In the Lund interpretation of 
fragmentation where the particles stem from the energy of the force field it is 
tacitly assumed that modifications of the perturbative state below and up to 
the scale of the hadronic masses should have no effects. We will find that it is 
necessary to take into account the coherence properties of the radiation in any 
modification. 

Although the methods presented in this paper are applicable to any multi- 
gluon state, for definiteness we will concentrate on the states obtained in e + e~- 
annihilation processes where we expect an original colour singlet (qq)-st&te to 
form and start to go apart, producing a set of bremsstrahlung gluons inside 
an essentially point-like region. We will also be satisfied to treat a single kind 
of hadron with mass m. Finally, in this paper we will not introduce gaussian 
transverse momentum fluctuations (which we expect in a tunneling scenario 
H) in the fragmentation process. We will investigate the influence of such 
fluctuations in a future publication. 

In Section ||, we provide a set of necessary formulae from the (1 + 1)- 
dimensional Lund Model. In Section ||, we consider the motion of strings con- 
taining internal excitations. We also provide a description of the coherence 
properties of QCD and some necessary formulae to understand the A'-curve. In 
Section |J, we define the properties of the string breakup in general. We find 
that the most general implementation of the Area Law leads to a highly in- 
tractable process and therefore in Section || we define another approach which 
mends all the problems with the earlier one. In particular, in Section 5.3 we 
show a close relationship between a differential version of this process and the 
A"-curve. In Section ^, we present a set of results from our method, followed by 
some concluding remarks on future work in section 7. 



3 



2 Some results from the (l+l)-dimensional model 



2.1 The Area Law 

The Lund Model contains a non-trivial interpretation of the QCD force-field in 
terms of the massless relativistic string with the quarks (q) and the anti-quarks 
(q) at the endpoints and the gluons (g) as internal excitations on the string 
field. It is assumed that the force field can break up into smaller parts in the 
fragmentation process by the production of new (gg)-states (i.e. new endpoints). 
A q from one such breakup point ("vertex"), together with a q from an adjacent 
vertex and the field between them, can form a hadron on mass shell. 

For the simple case when there are no gluons, the string field only corre- 
sponds to a constant force field (with a phcnomcnological size k ~ IGeV / fm) 
spanned between the original gg-pair. In a semi-classical picture, conservation 
of energy-momentum allows the creation of a new massless pair at some point 
along the field. The pair will then go apart along opposite light-cones, using up 
the energy in the field in between (in this way the confined fields will always 
end on the charges). In order that the hadron produced between two adjacent 
vertices should have a positive squared mass, it is necessary that the vertices are 
placed in a space-like manner with respect to each other. Consequently, time- 
ordering will be a frame-dependent statement (in any Lorentz frame the slowest 
particles will be the first to be produced, thereby fulfilling the requirements in 
a Landau-Pomeranchuk formation time scenario). It is possible to order the 
production process instead, along the light-cones and introduce the notion of 
rank so that the first rank hadron along the original g-light-cone will contain 
that q together with a q from the first vertex along the light-cone, the second 
rank hadron a q from the first vertex and a q from the next etc., cf. Fig|l]. Rank 
ordering is a frame independent procedure. It is of course possible to introduce 
a rank-ordering also from the end containing the original q. 

One obtains , §, @, [§, the unique process described by Eq. ([I]) from these 
observations and an assumption that the breakup process obeys a saturation 
assumption, i.e. that after many steps when we are far from the endpoints 
the proper times of the vertices will be distributed according to an energy- 
independent distribution. 

A particular feature is that if a particle with energy-momentum p = (p + ,p^) 
and with squared mass m 2 = p 2 = p+p- is produced between the two vertices 
with x = (x + ,x-) and x' — (x' + ,x'_) then we have (cf. Fig.^): 



P+ = K ( x + ~ = q+ - Q+ 

P- = k(x'_ — X- ) = q- — q'_ (2) 

Thus we find that on a flat string surface the difference between the vertex 
points will fulfil : 

(a: - x'f = -m 2 /K 2 (3) 



4 



Figure 2: Two adjacent vertices of the string breakup process. 
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Figure 3: A ladder diagram describing multi-particle production. 



Eq. (j|) implies that the (1 + l)-dimensional Lund Fragmentation Model may 
also be described by means of a multi-peripheral chain diagram as in Fig.|| or 
in an energy-momentum space picture as in Fig.^j. 

This is used in (5) in order to subdivide the Area Law process into steps 
in between the vertices. The energy-momentum conserving (5-distribution in 
Eq. (0) can be "solved" by introducing the momentum transfers {qj} instead 
of the hadron momenta {pj}- Then the mass-shell condition means that the 
hyperbolic angle between the vertices and the size of the area slit, exhibited in 



„n2 



Fig.|| are fixed by the squared sizes q 2 — —T, (q') 2 = —V and (q — q' 
The result is that Eq. ([[]) can be rewritten as a product of steps between the 



dP n ({ Pj },P tot )=l[K(T 



3^3- 



l, m 2 )dTj 
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Figure 4: The energy- momentum space picture (qj) and the spacetime picture 
(xj) respectively. 



, 9N exp(-%J\(r,T',-m 2 )) 

K(r,r',m 2 ) = jv y , 2V v ; ' - 
VA(r,r',-m2) 

A (a, 6, c) = a 2 + 6 2 + c 2 - 2afe - 2ac - 2bc (4) 

It is a remarkable fact that the transfer operators K can be diagonalised in terms 
of the eigenfunctions of the harmonic oscillator (those which are boost-invariant 
in a (1 + l)-dimensional space-like Minkowski space, in a two-dimensional Eu- 
clidean space they correspond to a vanishing angular momentum) <?„(r) with 
the eigenvalues solely determined by the squared mass of the hadrons produced 
in between: 

oo 

k{t, r', m 2 ) = ]T <? n (r)A„(m 2 ) 3 „(r') (5) 

71=0 

The eigenvalues A„ are analytic continuations of the harmonic oscillator eigen- 
functions to time- like values of the argument, ||. Useful representations of K 
and the eigenvalues A„ are : 

1 dz , b ~~ 2 - 2 



(6) 



^(r,r',m 2 )= / _ex P (--(zr + — ))s(V - (1 - z)(r + — )) 

J z 2 z z 

, / ?\ „r ,bm 2 . f 1 dz . .„ , 6m 2 , 

A„(m 2 ) = iVexp( — ) / — (l-z) n exp( ) 

l Jq z z 



We have introduced the positive light-cone fraction of the produced hadron 
z defined by (x + — x' + ) = zx + . It is straight forward algebra to prove that 
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Figure 5: The figure shows the two vertices and the regions I, II, III and IV 
described in the text. 

the shaded area, exhibited in FigJ| is given by the exponent \{zT + 2^) in the 
representation of the kernel K. We also note that the area obtained by summing 
the areas of the regions marked I, II and III in Fig.[| is equal to twice this 
area and that m 2 / z equals the sum of the areas of regions I and 77. There is a 
simple relationship between the two adjacent values of T in the representation 
of K: 

r' = (i-z)(r + ^) (?) 

Finally we note the identity (in easily understood notation) : 

(/)(///) = (II)(IV) (8) 

Eq. (§), just as Eq. (||), is only valid for a flat string surface and a particular 
consequence is that the region IV is non- vanishing for the (1 + l)-dimensional 
model. As the area of the region IV is proportional to (1 — z), the variable z 
must always be smaller than unity, i.e. there is a built in requirement that a 
typical step of the process can never use up all the available light-cone energy- 
momentum. 

Such a requirement also comes out of the following argument. Suppose that 
we would integrate dP n in Eq (^) over all possible energy momenta and then 
sum over all multiplicities. Due to Lorentz invariance, we will obtain a function 
R(s) which can only depend upon the total squared energy-momentum s = P 2 ot - 
If we pick out the dependence on the first particle and sum and integrate over 
all the rest we obtain an integral equation for the function R : 

R(s) = (B.T)+ / N— exp(— 6 — )R(s') 
Jo z z 

S ' = (l-z)(s-—) (9) 
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where {B.T) stands for "boundary condition term" and the variable s' is equal 
to the squared mass of all the remaining particles if the first hadron takes the 
light-cone fraction z (we note the similarity to the Eq. ([?])). The integral 
equation (|^) has an asymptotic solution R oc s a (with the parameter a being a 
function of N and bm 2 ,cf. Q) with the requirement : 

i d 7 m 2 

/ N— (l- z ) a exp(-6— ) = 1 (10) 
Jo z z 

Consequently while the exclusive formula for the production of a particular 
hadron with the light-cone fraction z is given by the Area Law, the inclusive 
probability to produce this hadron (irrespective of what comes after it in the 
process) must be weighted with R(s') / 'R(s) ~ (1 — z) a . Therefore, the well- 
known Lund fragmentation formula is given by the integrand in Eq. ( |Io|) and 
there is a power suppression for large values of the fragmentation variable. 

The formulae presented above correspond to an ordering along the positive 
light-cone, i.e. the variable z is defined as the positive light-cone momentum 
fraction of the particle. It is possible to redefine everything in terms of an order- 
ing along the negative light-cone, i.e. to introduce the corresponding negative 
light-cone component £ by writing (x'_ — X-) = C,x'_. It is straightforward to 
prove that 

9 2 

£ = 2 i r and 2 = 2TTF1 ( n 
m + zl m z + £1 ' 

and from this we find that the integrand in the representation of the kernel K 
can be reformulated from (z, V) — > (£, V) to exhibit the symmetry between the 
descriptions along the positive and the negative light-cone directions 

^6(T> - (1 - ,)CT + ^)) exp(-^zr + ^)) - 

^S(T (1 - C )(r' + ^)) exp(-^(Cr' + ^)) (12) 



3 The description of a multi-gluon string state 

The dynamics of the massless relativistic string is based upon the principle that 
the surface spanned by the string during its motion is a minimal surface. This 
means that the surface is completely determined by its boundary. In the Lund 
Model the string is used as a model for the confined colour force field in QCD 
and the above property then has the further important implication that the 
dynamics will be infrared stable, i.e. all predictable features from the decay of 
the force field should be stable against minor deformations of the boundary. 

For an open string a single wave moves across the spacetime surface and 
bounces at the endpoints. The wave motion is determined by a (four-)vector- 
valued shape function, which we will call the directrix, A. Thus a point on the 
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string, parametrized by the amount of energy a, between the point and (for 
dcfiniteness) the g-endpoint is, at the time t, at the position 

x(a,t) = \(A(t+-)+A(t--)) (13) 

We will from now on put the string constant n equal to unity in order to simplify 
the formulae. 

While the tension T = dx/da is directed along the string, the velocity 
v = dx/dt is directed transversely so that T ■ v — 0. The definition of a also 
implies that T 2 + v 2 = 1 (all the three-vector relations are valid in the local 
rest- frame). Together this means that the directrix function has an everywhere 
light-like tangent 



(14) 



The tension must vanish at the endpoints (a = and a = E to t) and this implies 
that the directrix must be a periodic function with the property 

A(£ + 2E tot ) = A(0 + 2P tot (15) 

where Ptot (Etot) is the total energy-momentum (energy) of the state. While 
according to Eq. (fl3|), the directrix A(t) describes the motion of the g-end, from 
Eq. ( |l5| ) it is evident that 

A q (t) = A(t + E tot ) - Ptot (16) 

will describe the motion of the q-end. Finally, if the string starts out from a 
point (at the time t = 0) then we must have the symmetry : 

A(0 = -A(-() (17) 

Using the Lund interpretation of the gluons as internal excitations on the string 
it is easy to construct the first half period of the directrix curve: it starts with 
the quark energy-momentum k\ and then the gluon energy- momenta {kj} are 
laid out in colour order and it ends with the q energy-momentum k n , for a string 
with (n— 2) gluons. In this way the g-endpoint will be acted upon by the colour- 
ordered excitations as they arrive in turn. From Eqs. ( |i~5| ) and ( |l~7j ) it follows 
that we obtain the directrix of the second half period by reversing the order, 
starting with the q energy-momentum and ending with the q energy-momentum 
(besides the translation this is the way the g-endpoint will move according to 
Eq. ©). 

The energy-momentum content in the string at a certain time t, between the 
point a and the g-end is given by : 

da'^ = \(A(t + a)-A(t-a)) (18) 
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3.1 The coherence conditions in QCD 

The properties of the directrices which are described above are common to all 
states of the massless relativistic string. On the other hand the use of the string 
as a model for the force fields of QCD may single out a particular class of all 
possible states. We will briefly discuss the conditions which correspond to the 
coherence conditions of bremsstrahlung radiation in a gauge field theory. 

Multi-gluon radiation is in general described by means of perturbative cas- 
cade models. In order to consider the properties of the states, we will make 
use of the ideas behind the Lund Dipole Model. To see the emergence of a 
multi-gluon state we start out with the following two basic results from QCD 
bremsstrahlung, (cf. §,§,§)• 

Bl The original (gg)-state emits bremsstrahlung radiation according to the 
ordinary dipole formula, i.e. there is an inclusive density of gluon quanta 
which is flat in rapidity y and the logarithm of the squared transverse mo- 
mentum k — log(fcj_) with a density given by the coupling, a = Ca s /27r: 

dn = adndy (19) 



B2 After the first ("hardest") gluon (gi) is emitted there are two dipoles 
available, one between the (qgi) and one between the {giq) with squared 
masses S12 = (k q + k gi ) 2 and S23 = {k gi + kg) 2 . The second gluon may 
be emitted independently inside the angular (rapidity) range of the first 
or the second dipole. We note that the two dipoles (each spanned by two 
light-like parton energy momenta) move apart in the total cms. 

After the emission of the second gluon, according to the Lund Dipole Model 
there are three independent dipoles for further emission, and so on. The masses 
of the dipoles will quickly decrease and therefore also the transverse momentum 
size of the emitted gluons. The cascades are stopped by some "virtuality" cutoff, 
either in the dipole mass or in the transverse momentum. The Dipole Model is 
implemented in the Monte Carlo simulation program ARIADNE || . 

In terms of the directrix, the original state is described by the light-like 
energy momenta of the original (qq)-j)ah:. After the emission of the first gluon 
the state, in space (in the cms), will be described by a (connected) triangle, 
cf. Fig.^, where the vectors k q = k\ and the vectors k g = k,2 constitute one 
of the dipoles and the ki and k q = k^ the second dipole, as discussed above. 
(We note that for light-like vectors the space part length is equal to the time 
component) Emission of the second gluon will then occur in between the two 
vectors describing the dipole and this evidently means that the second gluon 
vector will "cut off" one of the triangular corners, making the directrix into a 
quadrangle etc. (Note that already with the emission of two gluons, the vectors 
shown in Fig.^| no longer need to be in a plane) 

In general, there will be an emission of a set of "hard" gluons which will 
determine the general shape of the directrix. The remaining emissions will then 
make the directrix smoother and smoother as each new emission will correspond 
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Figure 6: The emission of the first and the second gluon respectively. 



to a gluon vector which cuts off an earlier corner. This means that the angle 
between the energy-momentum vectors of colour-connected partons becomes 
smaller (along the main directions, determined by the hard emissions) the longer 
the cascade continues. This is the way that the coherence properties of the 
bremsstrahlung radiation work, and it is sometimes referred to as the "strong 
angular condition" . Later on, we will find that this angular condition will play 
an important role when we consider the deviations in the partonic states which 
are allowed by the fragmentation process. 

4 The General Breakup of a String Field 

We will now consider the partitioning of a general string state at a point x(a, t) 
and after that define the most general process possible for the Area Law. 

According to Eq. ( |l3| ) the breakup will occur at the middle point \(A(£,\) + 
^(£2)) between two positions on the directrix, determined by £1 > £2 with 
t = + £ 2 ) and a = |(£i — £ 2 ), cf. Fig.^. There will be two parts of the 
string left over and we will now describe their motion after the breakup. 

The first part (to be denoted "(12)") can be described by the "new" directrix 
A±2 (£) = -4(C) f° r £2 < £ < £i- It will contain the energy-momentum P\2 — 
(Afa) - A(&))/2, cf. Fig.0, so that it can be continued as in Eq. @ with 
Ptot — > P\2- Starting it out from the position ^12(^2) = -Mfci), it is the "new" 
orbit of the (original) g-particle in the string (12). To obtain the corresponding 
orbit for the "new" <?-end we use the Eq. ( |l6|) to obtain A\2q = A\2(t + (£1 — 
^2)/2) — P\2- If we start at the "breakup" time t = ^(£1 + £2), it will evidently 
behave just as if we had adjoined the directrix A12 to the point cc((7, t) (note 
that ^i2(£i) = 2P 12 +^ 12 (£ 2 )). 

The second part (denoted "(21)") can be described as the remainder, i.e. 
^4 2 i = -A(£) with £1 < £ < 2E to t + £2 (noting that the original directrix is 
continued according to Eq. (|15|)) . The directrix .A21 contains the energy 
P21 = Ptot — P\i and is continued accordingly. To find the orbit of the produced 
q-end in the part (21) we adjoin the directrix _A 2 i to the point x(a, t) starting 
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Figure 7: The system (12) is shown as a string piece which is moving away, 
while the system (21) is translated to x\2- The second break-up point X34 is the 
middle point between A(^) and ^(£4). In this way, the particle produced in 
between X12 and X34 can be taken either as the second particle in the break-up 
of the original system or the first in the system (21). 



it at £ = £1. The (original) q will again move according to Eq. ([l6|). 

The most noticeable property is that for a multi-gluon force field the direc- 
trices of the two new string parts will not fulfil Eq. (p"7j). This implies that 
the endpoints of the new strings never meet but instead turn around each other 
so that each of the two states will contain angular momentum. It is only in 
the (1 + 1)- dimensional case that the two new string parts will have the same 
properties as the original string (though they will be scaled down in size and 
have different rest-frames due to the momentum transfer at the breakup). The 
conclusion is that the final state in a multi-gluon string fragmentation in general 
depends upon the ordering of the different breakups. 

4.1 A general process based upon the Area Law 

We may nevertheless devise a string breaking process on a "frozen" string sur- 
face, i.e. under the assumption that the string surface is once and for all a given 
device (this is also the way Sjostrand treated the problem exhibited above). 
Assume that the system (12) is a hadron on the mass-shell, produced as the 
"first" hadron along the original g-direction (i.e. the point A(&) is along the 
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Figure 8: The break-up points and the corresponding right and left-moving 
grains in the (a, t)-plane, as described in the text. 

q energy-momentum vector and A(£i) on the other side of the corner in the 
directrix between the q and the first gluon). Then we may choose two new 
points £3 > £1 and £4 < £ 2 to obtain the second rank hadron (34) with energy- 
momentum P34 = |(-4(^3) — -4(d) + -M,^) — -M^i)) with the production vertex 
at X34 = I (-4(^3 ) + Afa)). We may evidently continue this process across the 
string surface and we note that at every step it is only necessary to choose two 
numbers: (5£)+ — $3 — £1 and = £2 — £4. They must be chosen so that 

the corresponding J334 is the energy-momentum of a particle on the mass-shell. 
It is necessary to have a second condition, however, and we may then use the 
Area Law. 

For the (1 + l)-dimensional model the numbers (<$£)± evidently correspond 
to the light-cone coordinates z and £ in Section |2| The indices ± can be given 
a further meaning for the general directrix. According to Eq. ( |l3| ) there is a 
left-moving and a right-moving wave spanning the string surface and according 
to Eq. Jl5|) the fixed phase parts (from now on we will call them "grains") 
"stream" across and bounce back (thereby turning from left- to right-movers) 
at the q and q ends. 

We may describe the situation in a (a, <)-plane, cf. FigJ|, with the right-and 
left moving grains moving along fixed "light-cone lines". The whole string is, 
at the time t — 0, gathered into a single point at the origin. The boundary 
values at a = and a — E to t are the directrix A(t) and the orbit of the q, i.e. 
Aq(t), respectively. A breakup point x = x\i is reached by the meeting of the 
left-moving grain indexed £1 and the "right-mover" £2 and the breakup at £34 
by the left-mover £3 and right-mover £4. The energy- momentum P34, streaming 
into the hadron produced between x\2 and X34, is evidently i(^l(£ 3 ) — -4(£i)) 
from the left and ^(A(£,2) — ^(£4)) from the right. If we use ± as indices for 
the left- and the right-movers respectively we obtain a useful parametrisation. 
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With regard to surface areas there is an obvious (scalar) surface element 
defined by the surface spanned by the (incremental) grains coming from the left 
and the right. 

dA oc dA+ ■ dA- (20) 

Using this we may easily calculate the area "below the vertices" a; 12 and £34 (cf. 
Fig.|J) as they are defined above. They are the correspondence to the squared 
proper times T considered in Section ^ (in the Fig. |^ they would correspond to 
the areas I + IV and 77/ + IV). 

Ti2 = x\ 2 and T 34 = x\ A (21) 

So, for the vertices shown in Fig|3, these areas are just the squared proper times 
of the vertices. But we note that they can not be expressed as the products of 
light-cone components or even in terms of products of sets of left- and right- 
movers solely, for general multigluon strings. 

The area below the vertices (if we neglect the squared mass term in the Area 
Law, which can evidently be included in the normalisation constants N in Eq. 
(0)) can be expressed as 

A1234 = x\ 2 + £34 - x\ A (22) 
(this is the correspondence to the area called / + III + IV in connection with 



Eq. (|8|)). We note that the area corresponding to IV, i.e. x\ A in Eq. (22), may 
vanish in this case, i.e. the left-mover indexed 1 and the right-mover indexed 4 
may meet "before" t — 0. This means that the condition in Eq (||) can not be 
met for the general situation. 

We may nevertheless include these considerations to devise an exact version 
of the Area Law in Eq. (|l|). We write, in an obvious way, the differentials 
d 2 p — > (dp+ ■ dp-) and include the mass-shell condition by a ^-distribution as 
in Eq. (|l|). The result is, however, very complex to treat and for our purposes 
there are three shortcomings: 

1 If we treat the process in an iterative way, at every step, just as we men- 
tioned above, there are two numbers (#£)± to be solved for. One is fixed 
by the mass-shell condition and the other can be fixed by the area expo- 
nential; but the distribution functions are very complex (cf. the remarks 
under 2 below) . In principle everything can be done by Monte Carlo sim- 
ulation methods but even with very fast computers and the best possible 
computer routines this model will always be much slower than a process 
where there is only one number to be solved for, at every step. 

2 Due to the complexity of the formulae there is no way we can obtain 
useful analytical tools to study the behaviour of the correspondences to 
the transition operators as in e.g. Eq. (||). This also implies that there is 
no way that we can define the total integral and sum over the produced 
hadronic states (as we could do for the (1 + l)-dimensional model in Eqs. 
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(JsJ) and (|To|) ) so as to be able to obtain the correct weighting in every step 
of an inclusive cascade. 

3 Although this implementation of the Area Law is an "exact" procedure, 
fulfilling at every step the mass-shell condition and taking the "true" area 
below the vertices into account, it nevertheless starts out from the as- 
sumption that the partonic state is so well defined that there are reasons 
to implement an "exact" theoretical fragmentation procedure on it. This 
is in general not the case, i.e. the partonic state defined by perturba- 
tion theory is only well-defined down to some virtuality cut-off as we have 
discussed before. 

We use quotation-marks on the word "exact" in order to indicate that the use 
of a "frozen surface" is in itself an assumption introduced to mend the problems 
obtained from the use of classical physics and neglect of the subsequent angular 
momentum production. In the third item above, we are coming back to the 
statement at the beginning of the section, that infrared stability should imply 
that minor deformations of the boundary (in this case changes in the partonic 
state as defined by the directrix A up to the order of the hadronic mass scales) 
should not be noticable in the results. In the next section we will define a 
procedure which will mend all the three shortcomings mentioned above. 

5 The Lund String Fragmentation as a process 
along the directrix 

In this section, we will consider a fragmentation process which is defined along 
the directrix and does not suffer from the three problems considered at the 



end of Section 4.1. In order to exhibit the idea we will start with the (1 + 1)- 
dimensional model and rewrite it in a useful way. After that we will extend it 
to the general case. In particular, we will allow modifications of the order of 
the hadronic mass-scale in the directrix in accordance with the discussion under 
point 3 above. 

5.1 The directrix process for the (1 + l)-dimensional case 

In this case the directrix contains only two directions, given by the q energy- 
momentu m ve ctor (to be called A+ in accordance with the notation introduced 



in Section 4.1 ) and the q energy-momentum (-4-). A vertex point Xj, obtained 
after the production of j hadrons from the g-side, pi,...,pj is then described 
(with respect to the origin) by 

hA+i+A-i) (23) 



2 



We also know that 



j 



1 (A +j -A. j ) (24) 



i 



2 
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Figure 9: A description of the fragmentation process along the directrix with 
notations according to the text. 



Using the symmetry of a directrix passing through a sing le point (Eq. @) we 
may find another point on the directrix with the property 



A- 



(25) 



We will from now on drop the indices ± on A+ and B~ but we note that 
they do describe points on the same directrix in accordance with the left- and 
right-mover notation in Section 4.1. While A- goes "backward" for increasing 
j-values, B follows the g-direction. 

We may now consider the hadron energy momenta to define a curve from 
the origin "along the directrix" such that after j steps it has reached the point 



(26) 



while the difference between the point Aj on the directrix and Xj is given by Xj 
in Eq. (^3|). The production of a new particle Pj+i then corresponds to choosing 
two new points X/+i and Aj+i (along the directrix) such that, cf. Fig.^ 



Xj + i — Xj = Pj+\ and Aj+i — Aj = fcj+i 



(27) 



We also obtain a new "vertex" vector a^+i by the identity ( cf. Fig.g): 

Pj +X + Xj + i = Xj + kj+i (28) 

In this way the vertex vector fulfils x^+i = Aj+\ — -Xj+i just as Xj = Aj — Xj. 
We have then arranged it so that the hadrons are produced along a curve, the 
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X-curve, from the origin and the vertex vectors are the connectors for this curve 
going from the produced particle to the directrix. Before we consider the Area 
Law in this situation we note the symmetry between a reversed process and the 
process described above, i.e. when we go from Xj to Xj + i producing pj+i by 
the use of a part kj+i of the directrix along A. 

To see the reverse process we note that the vector Xj can just as well be 
reached by taking the difference between the point Xj on the hadron curve, Eq. 
(p6|), and "the backward point" on the directrix Bj. 

Xj = i {Aj -Bj) = Aj - X 3 = X 3 - Bj (29) 

Using this we could evidently consider the production of the particle Pj+i as a 
step from Bj to By+i = Bj + lj (cf. Eq. (|27|)) such that we have in correspon- 
dence to Eq. @ and Fig.| 

Vj+i + Xj = Xj+i + lj+i (30) 

In order to formalise the determination of the particle energy momentum p, we 
may then in "the fc-process" (along .4) assume that we know the starting vertex 
vector x, connected to a point Ap. We may then choose a piece k along A (of 
a size to be determined) and then define the other light-cone direction in the 
plane determined by (x, k) by 

x 2 

l = x -k— (31) 
2xk v ; 

The vector p will be described in terms of (fc, £) as 

p = z £+ = zx +-(l-— - (32) 
2 2 xk 

with the requirement that the particle should be on the mass-shell 



..i 2 

p 2 = m 2 = zkx i.e. kx — — (33) 



9 2 ; • i m 

•' = zkx i.e. kx = ■ 

z 

From Eq. (p8|) we obtain the new vertex vector x' by 



^ = (1 -*)*+!(! + ££) 

2 

{x'f = (1 - z)(x 2 + xk) = (1 - z)(x 2 + — ) (34) 

In Figjl^ we show both the production as described in Section || and the fc- 
process described above. In particular we note the two areas exhibited. It is 
straightforward to prove that 

a = A r + r' - r (35) 
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Figure 10: The figure shows both the production as described in Section g and 
the ^-process described in the text. 

with the obvious definition of e.g. T = x 2 . The Area Law in Eq. ([[]) will 
be fulfilled both by the use of the area in Fig.[l] and the area in between the 
directrix and the hadronic curve (the X-curve) if we choose the variable z from 
the fragmentation function in Eq. ( |To| ) and apply it as in Eqs. (^2|)-(^). This 
is so because the difference between the areas A and A r for the single step 
exhibited above will vanish when we consider the whole process. 

It is also obvious that we may define an ^-process similar to the /c-process 
we have discussed above. We just write zi = M and introduce the variable ( 
such that ife = £fc with £ chosen such that m 2 /( = tx' . Actually we obtain 
the same process (although "in the opposite order") under the assumption that 
we start at x' and chose I along the 2?-part of the directrix with the variable ( 
in accordance with Eq. (|ll|). In this way the "backward" variable £ evidently 
obeys the same distribution as the "forward" variable z and the Area Law is 
fulfilled. 

We have reached a situation where a particle production step starts from a 
knowledge of a vector x connected to a light-cone-direction. Then we choose a 
light-like vector k such that Eq. ([53|) is fulfilled with a z-value stochastically 
chosen from the fragmentation function in Eq. ([To]) . After that we construct the 
particle energy-momentum and a new vector x' according to Eqs. ( |32| ) and (p4|). 
We may start out choosing the "first" x- vector (xo) equal to the q (light-cone) 
energy-momentum and then the (1 + l)-dimensional model is defined. 

5.2 The directrix process around a gluon corner 

For a multi-gluon directrix however, we will reach situations when there are "cor- 
ners" , i.e. where the directrix changes direction between two colour-connected 
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partons. Thus besides the vertex vector x, there will be a remaining light-cone 
vector (to the corner) c and a chosen number z such that 



9 

m 



cx < (36) 



z 



We will now consider ways to pass over such a corner. Assuming that the 
directrix continues on the other side of the corner we may choose a light-cone 



vector k which fulfils Eq. ( |33| ) (we note that according to Eq. (36) k can not 
be chosen along c) and another light-cone vector d such that 

k + c' = c + Q = Q (37) 

with Q, a vector pointing from the corner to a new point on the directrix, and 
Q describing that point from the point on the directrix where the vertex vector 
x is connected. 

We note that with this definition the suggested way to pass the corner cor- 
responds to an exchange of one "dipole" , c + Q in the process for another k + d, 
i.e. to make a detour onto a "new" directrix which is "close" to the original 
one. The choice of the new dipole corresponds to a definition of what we will 
mean by the allowed "minor deviations" of the original perturbative state. 

In order to investigate what can be allowed under these conditions we will 
exhibit a few possible choices to fulfil Eq. (^7|). We will then find that all 
possible choices are not allowed if we would like to obtain reasonable results for 
the final state hadron distributions. A first set of possibilities is to choose d as 
a new light-like connector that will bring us "closer" (than the vector c) to the 
original directrix in the next production step. One may then hope that after 
another step in the process the new connector becomes even "smaller" so that 
we are back at the original directrix after a few steps. 

CI The vector k can be chosen in the plane spanned by the vertex vector x 
and the directrix Q (then both d and k are fully determined). 

C2 A more general choice is to take the vector k in the three-space spanned 
by x, Q and the remainder vector c. In this case it is necessary to provide 
a second condition to determine k and d. In the rest-frame of the vertex 
vector x the space parts of the vectors c and Q span a triangle with the 
vector Q as baseline, cf. Fig.[llJ (For the case exhibited we assume that 
also Q is light-like. The situation can be easily generalised to the situation 
when there are one or more "corners" also on Q.) The remaining degree 
of freedom will then be fixed by an angular condition. One such condition 

is indicated in the Fig.pl]. The two triangles (c, Q, Q) and (d, k, Q) in this 
construction are chosen congruent. This is a kind of minimal choice for 
the size of Q in this three-space. It is not difficult to see that in general 
the new connector will be smoothly connected to the original directrix. 
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Figure 11: A possible construction (C2) for passing a gluon corner. The figure 
shows the space parts of the vectors, in the rest-frame of the vertex vector x, as 
described in the text. 



For the case defined by the condition CI we obtain that k and d will be 
directed along the two light-cone directions in the plane spanned by (x, Q): 

2 2^{Qxf - Q 2 ^ 

(In Fig.|ll| they would be directed along the Q-direction with the lengths (Qo ± 
|Q|)/2, respectively). The vector Q contains one degree of freedom which can 
be fixed by multiplying with x in Eq. (|38|), and using Eq. d&jj). 

This means on the one hand that the situation corresponds to using up the 
smallest possible segment of the directrix in the construction according to Eq. 
(|37|). On the other hand we will obtain a very sharp "bend" on the new directrix 
segment, defined by the dipole (fee'). While the first result is desirable from an 
economy point of view, the second implies that we will break the coherence 



conditions for the directrix (the angular ordering) as discussed in Section 3^ . 

The fragmentation distributions obtained with this choice are not very en- 
couraging. After we "pass" the gluon corner the density of hadrons is much 
higher than the density before the corner. We find that even for a small gluonic 
excitation a break of coherence may result in the production of too many hadrons 
over a large rapidity range. 

One might hope that the situation may be mended by a suitable angular 
choice. This is so for the situation described in C2 for a large part of the 
events and for simple gluonic configurations. In Fig.|l8|we show a Monte Carlo 
simulation of the (inclusive) final state hadron rapidities in a state containing a 
single gluon placed at the centre with a transverse momentum of 15 GeV. Even 
though for most fragmentation events the choice C2 performs well, on the level 
of a few percent it does happen that the directions of the connecting c' and 
the new vertex x' are such that we obtain a situation close to the one for CI. 
The new directrix (possibly after a few further steps in the process) will all the 
time come closer to the original one; but at the same time we obtain a sharp 
bend farther and farther away from the original gluon corner with each step. 
This causes too many particles to be produced over a large rapidity range (cf. 
Fig®. 
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It is possible to make further changes in the procedure with more sophisti- 
cated angular choices, but we will instead go over to another and more successful 
choice, which we will call C3. 

C3 A natural choice for the vector d is to make it into the /c-vector for the 
next particle production in the process, i.e. chose a new stochastic value 
z' and put d — k' with 

k'x' = — (39) 

z 

in terms of the new vertex vector x' determined from (z,k,x) according 
to Eq. (|34|). In this way we are evidently getting back to the original 
directrix as fast as possible. 

Wc may define the vector k still in the space spanned by (c, x, Q): 

k = ac + 0x + 7Q (40) 

We obtain a solution for the coefficients (a, (3, 7) from the requirements 

k 2 = 

777 ^ 1 

kx = — and kQ = -Q 2 (41) 

(these conditions also imply that k' is light-like), and an equation for the vector 
Q along the directrix using Eqs. @ and 

^ = k'x' = (Q - fc)((l - Z )x+ k -(l + ^)) (42) 

In the fragmentation distribution there is a factor (1 — z) a (stemming from the 
density of states according to the discussion around Eq. (10)). This implies 
that we may not use up all the energy-momentum along one of the light-cone 
directions in a single step. Small values of z correspond to C — 1 according to 
Eq. (11), i.e. to using up almost all the energy-momentum along the other light- 
cone direction in a step. Such values are suppressed by the area exponential in 
the distribution. Nevertheless it may happen that we obtain a very small z-value 
[z < 0.1), corresponding to a large (but not infinitely large) area. Sometimes 
for the general multi-gluon directrix however, there might be no possibility to 
accomodate such a large area due to energy-momentum conservation. Then the 
choice C3 can not be done. Instead of introducing a general projection of the 
two-particle distributions for all such cases, in the Monte Carlo simulation we 
have chosen to go back to the beginning and start a new fragmentation event. 
Such situations, when we can not use the choice C3, are infrequent enough to 
let this "restart strategy" be usable. 

In order to test the influence of this feature on the model we investigate the 
influence on the single particle and two-particle distributions below in Section 
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^ and we find them negligible. In order to show that our changes in the direc- 
trices actually are of the hadronic mass scale we have defined a local "bending 
parameter" closely related to the k± variable used in the ordering of the Lund 
Dipole Model, cf. Eq. @ 

2 _ 2(xfci)(fcifc 2 ) . . 

fc±2 - x(h + k 2 ) + k,k 2 

From the results in Section ^ we conclude that we are at most making local 
modifications on the order or below the hadronic scale in our directrices. 



5.3 A Differential Process and its Relationship to the Gen- 
eralised Rapidity 

There is a direct connection between a differential version of our hadronisation 
process and the <Y-curve that was referred to in the Introduction, In order 
to investigate this process we consider the limiting situation for a vanishing 
mass parameter. Then the distribution function will obviously develop a pole 
for z — > 0. We will assume that the model is defined by the incremental step 
size dz with the ratio m/dz — > tuq. The corresponding incremental k vector will 
be called dA and it will fulfil the mass-shell condition (we will use the notation 
qp instead of x for the vertex vector in connection with the differential process) 

o 

qpdA = — > dzm\ (44) 

From the model formulae for the change in qp and the particle energy-momentum 
p (Eqs. (|3^) and (|J)) we obtain the following differential equations defining a 
curve to be called the "P-curve 



dA _ q% 
2 1 m 2 
dA g| 
2 { ml 

Firstly we note that from the sum and differences of the Eqs. ((45[) we obtain 



dV = dzq P + — (1-^r) 



dq P = -dzq P +— (1 + ^) (45) 



V + qp =A 

V -qp=C (46) 
where the vector C has a light-like tangent just as the directrix A: 

dC = 2dzq P - dA% (47) 
m 2 

In this way the P-curve goes in between two curves with everywhere light-like 
tangents and the vector qp connects to both of them. 
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Figure 12: A string with a q(ki), q(k^) and a single gluon excitation (fe)- The 
figure shows the connected region around the gluon "tip" . 



The vector is time- like and quickly approaches the length mo- To see 



that we multiply the second line of Eq. (45) by qp and obtain 

dq P = dz(-qp + ml) (48) 
This means (remembering that dz = qpdA/rn^) that we may write: 

q 2 p=ml(l-T p 1 ) with T P = exp(/(^)) (49) 

J m 

We have then assumed that qp — and Tp — 1 at the start of the process. We 
note that the integrand in the exponent for Tp is an area, more precisely the 
area betweeen the directrix and the 'P-curve. This is very similar to the results 
we obtained in || for the generalised rapidity and we will now briefly connect 
to these results. 



5.4 The A"-curve and its Properties 

The distribution in Eq. (|l|) contains two terms, the phase space and the expo- 
nential area suppression. In order to obtain a large probability it is necessary 
for a given total energy-momentum on the one hand to make many particles on 
the other hand to make them in such a way that the area is small. The obvious 
compromise is that the decay region is around a typical hyperbola with an av- 
erage squared distance to the origin < T >= Tq. The length of the hyperbola 
is proportional to the available rapidity range for the final state particles, i.e. 
Ay = log(s/ro) with s the squared cms energy. 

In the Lund Model interpretation for a string with a single gluon excitation, 
there will be two parts of the string; one spanned between the q and the g and 
one between the g and the q. Each of them should break up in a similar way as 
the single string region described by Eq. (0). Besides that there will be a few 
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particles produced in the connected region around the gluon "tip" , cf. Fig.p2|. 
If the energy-momenta of the partons are kj, j — 1, 2, 3 (with indices 1 (3) for 
the q (g)), then there will be two hyperbolic angular ranges, (Ay)i2 and (Ay)23- 
The sum of these ranges is 

A = (Ay) 12 + (Ay) 23 = log(^) + log(|^) 

= log(f)+log(^) (50) 

r 4r s 

Here Sji = (kj + kg) 2 and s = S12 + S23 + S13 and the factors 2 are introduced 
because only half of the gluon energy-momentum goes into each string region. 

The quantity k\ = s\ 2 s 2 2,/ s is a convenient (and Lorentz invariant) approx- 
imation for the transverse momentum of the emitted gluon. From Eq. ( p0[ ) 
we conclude that after the emission of a single gluon the phase space is in- 
creased from the single hyperbola result above by an amount corresponding to 
a "sticking-out tip" of length given by the logarithm of the emitted transverse 
momentum. In conventional notions this corresponds to the "anomalous dimen- 
sions" of QCD, i.e. the emission of a gluon increases the region of colour flow 
inside which more gluons can be emitted and hadronisation can take place. The 
whole scenario is easily visualised and used in the Lund Dipole Model, [|J, and 
the corresponding Monte Carlo simulation program ARIADNE, |J. 

It is straightforward to see that if there are many gluons then there is a 
corresponding quantity, a generalised rapidity A ~ log(Jl Sjj+i) stemming from 
the hyperbolae spanned between the colour-connected gluons. But we note 
that this is not an infrared stable definition. We will now provide a convenient 
generalisation. 

A closer examination of the region around the tip of a gluon reveals that there 
is a correction corresponding to a connected hyperbola in the region (k\, £3) be- 
tween the "endpoint" of the hyperbola in the region spanned between (k\, k 2 /2) 
and the one spanned between (k 2 /2,k 3 ), cf. Fig.[t2| In formulae we obtain for 
the average hyperbolae 

(aifei + i/?ifc 2 ) 2 = To and ( 73 fc 3 + ^f3 3 k 2 ) 2 = T 

(a 2 fci + l2 h + \k 2 f = T (51) 

with the ranges 1 > a x > 2T Q /s 12 , 2T /s u > a 2 > 0, 2r /s 12 < ft < 1 
and similarly for the other variables. The length of the two hyperbolae in the 
segments (ki,k 2 /2) and (k 2 /2,k-z) are then given by Eq. ( pOQ but the third 
hyperbola provides an extra contribution (in the appropriate limit S13 ~ s) 
equal to log(l + 4r s/si 2 S23). Then the total (generalised) rapidity length 
becomes 

Ai23-log(-^ + S§) (52) 



T (2r 



This is evidently a nice interpolation between the situations with and without 
a gluon on the string and it is also an infrared stable definition of the notion 
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of rapidity length. Eq. (|52j) is noted in [|2| and led us to introduce a functional 
defined on a multigluon string directrix. 

We may firstly define a set of connected integrals, 

In = / ds 01 ds 12 ■ ■ ■ ds nE (53) 



with the easily understood notation (cf. Eq. ([18[) ) Sjj+i — (A(£j) — -4(£j+i)) , 
i.e. it is proportional to the squared mass between the points £j and along 
the directrix. By performing the integrals, we obtain that the argument in the 
logarithm in Eq. ( [52"|) is given by the sum Ii/To + ^/(ZTo) 2 and that we may 
in general define the functional T by : 

T=l + Y — l " (54) 

as a suitable generalisation for any string state. For a finite number of partons 
N the terms with n > N will all vanish and we also note that the highest degree 
term will always have the generic form : 

o S12 S23 s N - 1N 
Ami 4m o 4m o 

We also note that for a finite total energy E the contributions for very large 
degrees will become smaller and smaller compared to the scale mo- 

In order to study the functional T it is suitable to introduce a varying value 
£ instead of the total energy E in the connected integrals. It is then evident 
that the functional T(£) will fulfil the integral equation : 



ds(tiA'), 



We will also introduce the vector- valued function qr(0 together with T so that 
we have 

dA,(QT(e) 

lTu{t,) = J^Tj 

T( =l + f^»T(n (57) 
Jo m o 

By differentiation and integration we obtain the results 
T = exp( / g ) = ex P( A (0) 



o 



q ^) = mUl-T- 2 (0) (58) 



The generalised rapidity A corresponds to the result in Eq. (|5^) for the sim- 
ple case described above and it provides an infrared stable definition for any 
multigluon state. Further the vector qx is time-like and will quickly approach 
the finite length mo. 
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5.5 The Correspondence between the A"-curve and the V- 
curve 



The similarity between these results and the results obtained for the P-curve 
in Eqs. ( [48| ) and ( [49| ) are obvious (besides the power in T _1 ) and there is a 
correspondence to the Eqs. ( |45| ) and ( phf ) also. The interpretation of the X- 
curve (as worked out in cf. also ||) is that there is a vector valued function 
X^(X) conveniently labelled by A such that 

X + q T = A 
dX 

lx =qT 

dq T . dA 

~d\ = qT ~d\ (59) 
i.e. the vector qr is the tangent to the curve defined by X such that it reaches 
to the directrix. There is no direct correspondence to the £-curve (unless the 
vectors q has reached its asymptotic length mo). 

It is useful to calculate the results from the differential equations for the 
case when there is a finite length light-like vector kj in the directrix. By direct 
integration we find that if we have the vector q^j then "after" application of the 
parton energy- momentum kj we obtain the vector Qtj+i and will take a step 
along the X-cmve equal to 5Xj 

(l+7i), 
QTj+i = IjQTj + g fc J 

v = 77^7 ( 6 °) 

(we also note that the products of the jj is equal to T _1 ). The correspondence 
to these results for the "P-curve are 



qp = lip h 2 



SV = (l-j)qp + 



(1 + ^?) 

v mi f 

(i 



2 

SC = 1 = 2(1- j)q P -^-k 
1 

1 _l_ g-pfe 



and (Tp)- 1 ^!]^ ( 61 ) 



This means that while the length of a step along the X-curve depends upon q\ 
the corresponding step length along the P-curve is 



{5Vf = Mf = (1 7j) "°» (62) 

7j 
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Figure 13: Two vertices placed symmetrically (around the rapidity y=0) at 
mo(cosh(<%/2), sinh(5y/2)) and mo(cosh(5y /2), — sivh(5y/2)) respectively, with 
a distance 2mo sinh(i5y/2)). 



This has a very simple meaning for a step length Mj along a hyperbola with 
the parameter mo, cf. Fig.jH?]. If the vertices are placed symmetrically around 
the rapidity y = at the positions mo(cosh(Jy/2), ± sinh(<5y/2)) then the step 
length is evidently 2mo smh(5y/2)) which should be compared to Mj. Then if 
we square them and define jj = (1 — Zj) — exp(— Sy) the relationship in Eq. 
( |62[ ) will ensue. 

The identification of 7 and (1 — z) is actually a very general feature of the 
process (we will come back to this aspect in the future). At this point we only 
note that if the lengths (proper-times) of the two adjacent vertices are equal 
r = T' = ml =< T > then from Eq. (Q) we obtain 

n (1 — z)m 2 
m o = J (63) 

which evidently coincides with the result in Eq. ([32]). Therefore, for this par- 
ticular value of z, the fragmentation process will all the time proceed along the 
connected hyperbolae. It is in this way that the P-curve can be considered the 
average of the hadronic X-curve. 

In order to provide a precise relationship between the vectors qp and qx as 
well as the functionals Tp and T we make use of an interesting relationship for 
the A"-curve which is derived in ||. If we define the (1 + 4)-dimensional vector 
(Qp = TqT/j,/mo,T) (which has a length in the (1 + 4) dimensional Minkowski 
metric equal to Q 2 — T 2 — —1) then the differential equation for qT can be 
rewritten as 



dT = QdA and dQ = TdA (64) 

i.e. as a group of special rotations in this space (corresponding to a subgroup 
of SO(l,4)) which are defined by the incremental changes along the directrix 
curve. The corresponding relationship for the "P-curve is, from the Eqs. ( ^i|) 
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and @ 



d(qpTp) = dA(T P - -) and dT P = q P dAT P (65) 

so that the (l + 4)-dimensional vector (2qpTp, 2Tp— 1) has both the same length 
and fulfils the same differential equations with respect to incremental changes 
along the directrix as (Q,T). In we characterized the different directrices 
by means of the eigenvalues of the transfer matrix along the directrix and we 
found a simple and useful method to relate them to the generalised rapidity 
A = log(T). We will come back to this formalism in a future publication and 
show the significance of the extended space. 

In conclusion we have found a differential version of our stochastic process 
which corresponds to a curve of connected hyperbolae along the directrix func- 
tion. The area in between the curve and the directrix (scaled by the single 
parameter m^) has the interpretation of a generalised rapidity measure for the 
multigluon case. It is also related to a group of rotations with the incremental 
steps along the directrix as the generators. The curve may be interpreted as the 
average hadronic curve stemming from the hadronisation of the given directrix. 



5.6 The Reverse Problem, to Find the Directrix from the 
Hadronic Curve 

We will end this section by solving the reverse problem to the hadronisation 
process, i.e. to exhibit to what extent we can trace the directrix from a knowl- 
edge of the hadronic curve, which we will call the X-curve in accordance with 



the notation introduced in Section 5.1 



We will then assume that the X-curve is defined by the hadronic energy 
momenta {pj}, ordered and laid out according to rank. (Note that this can not 
be done for a real event obtained from an experiment, since in general, rank is 
not an observable.) We will concentrate on the production of the hadron pj, 
produced between the vertex vectors Xj-i and Xj with the directrix segment kj, 
and an appropriately distributed stochastic number Zj . According to Eq. (Eq) , 
in order to construct kj it is sufficient to know pj and the difference vector 

( Xj - Xj-i) = ejpj (66) 

It is straightforward to solve for pj in terms of pj and Xj-\ 

p j = (67) 



2™2 
3 X 3-1 



The sign ej should be positive or negative depending upon whether m 2 / Zj is 
larger or smaller than ZjX 2 _ 1 (it is useful to note that 2(pjXj-i) — (m 2 /zj + 
ZjX 2 ^-^)). Therefore if we prescribe the first vertex vector xq (this is always cho- 
sen in our process as the original q energy-momentum vector) then the directrix 
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Figure 14: The x 2 distribution of the production vertices in our method for 
different values of the hadron mass, together with a plot of the F distribution 
in Eq.@ with a=0.5 and b=0.7. 



vectors as well as the vertices are determined recursively up to a sign: 

k j = /'.. + ' iPj 
xj = XJ--L + ejPj (68) 

It is evident that the other sign will determine the corresponding lj. We note 
however, that neither the necessary rank ordering of the hadrons, nor the colour 
ordering of the directrix gluons, are experimental observables. Therefore in this 
form the result has solely a theoretical meaning. It is also necessary to take 
the possible transverse fluctuations in the fragmentation process (mentioned 
before as stemming from tunneling) into account before any observables can be 
presented. 



6 Results 

In the previous section, we have proposed a method of fragmenting a multigluon 
string based upon the Area Law. It can be described as an iterative stochastic 
process along the directrix curve, which fulfils the Area Law at every step. We 
have discussed a few different approaches for passing a gluon corner, cf section 



5. 4 and from among these we have picked a preferred one. In this section we 
present a few basic results from this preferred solution and we also compare it 
with JETSET and the other approaches we have discussed. 

We start with comparing the distribution of x 2 from our model with the 
distribution of the proper time of the break-up points derived inside the (1 + 1)- 
dimensional model S . 

H(T) = NT a exp(-br) (69) 
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The result is presented in Fig.|lJ. As mentioned in the discussion at the end 
of section 5.2, sometimes certain values of z can not be accomodated when we 
pass the corner according to choice C3. Our Monte Carlo simulation program 
handles this by starting a new event from the beginning whenever this occurs. 
One may then expect that this might cause changes in the z-distribution and 
also in the T-distribution that we obtain. This is, however, not the case as can be 
seen from Fig. (14) and Fig. (15). Further these distributions are as they should 
be according to the Lund Model formulae. The T-distribution is completely 
independent of the mass value used just as in the standard, 1 + 1-dimensional 
Lund Model. 
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1.4 
1.2 

1 

% 0.8 
0.6 
0.4 
0.2 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

z 

Figure 15: z values used in particle production in our method (solid line) 
compared to the Lund fragmentation function f(z) (dotted line). Though z value 
generation is according to the distribution f(z), not every value produced is used. 
This plot shows however, that the rejection of z values in the present method is 
unbiased. 




To examine further effects in connection with our "restart strategy" , we de- 
fine the following measure Ay, of the correlation between two adjacent particles 
in the production process 

It corresponds to the difference in rapidity between the particles in the (1 + 1)- 
dimensional model. This distribution would be altered if we rejected particular 
correlations between successive values of z (e.x. a large value followed by a small 
value etc) more often than others. 
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Ay 

Figure 16: The Ay distribution, defined in eq.([70|), from our method (solid 
line) compared with the distribution obtained using only the Lund symmetric 
fragmentation function (dashed line). 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



k t2 2 

Figure 17: Plot of d ^ > vs k\ 2 measured in (GeV) 2 . Modifications of the 
directrix are on the scale of hadronic mass. 
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In Fig. 16 we compare our result with the result one obtains using only the 
Lund symmetric fragmentation function. There is no noticeable difference be- 
tween the two curves, suggesting that the restart procedure (involving possible 
rejection of some z values) is effectively insensitive to such correlations between 
adjacent values. 

In Fig.[f7| we show the distribution of our "local bending" parameter, A:^ 2 , 
defined in eq.(^). This figure illustrates that the modifications of the directrix 
that we make when we apply our method are of the order of the hadronic mass 
scale. For values of k\ 2 > m 2 there is a fall off faster than a gaussian. 




Figure 18: The figure to the left shows rapidity distributions for different choices 
for the vector k. The figure to the right shows rapidity distributions for a 3 jet 
event, obtained with the preferred method (C3), for different gluon energies. 
Note that the y axis is plotted on a log scale here. 

In the left-hand side plot of Fig.|l8| we summarize the rapidity distributions 
obtained using the different approaches for passing a gluon corner, C2 and C3 



described in section 5.2, and compare them to a curve produced using PYTHIA 
in one simple situation: a perturbative string state consisting of a quark, an 
antiquark and a single gluon. We conclude that modifying the directrix leads 
to the production of too many particles over a large range of rapidity if the 
modifications are not " smooth" . In the right-hand side of Fig.[l8| we show a few 
rapidity distributions obtained from using our preferred method for different 
energies of the gluon. 
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Figure 19: The rapidity distribution (left) and mean |p^ adrons | (GeV) as a 
function of rapidity (right), for a string state of A measure 70 (which is above 
the average A for the energy and p™' used for this plot: y/s — 2000 GeV , 
p™< — 2GeV ) generated by ARIADNE, fragmented using our method (solid 
curve) and PYTHIA (dashed curve). 



In Fig. (|19|) we show an example of inclusive rapidity distributions obtained 
by JETSET/PYTHIA and our method. For the comparison, we have used an 
arbitrarily chosen partonic event taken from ARIADNE. A close examination 
will show that in general there is a small multiplicity difference in the gluon 
jets. We have traced that to the fact that the two to three particles with the 
largest energy-momenta in the jets are faster according to our method compared 
to JETSET. It is impossible to compare to experimental data because firstly 
in this case we have only made use of a single partonic event. And secondly, 
at present we have not introduced the "fragmentation transverse momenta" . 
Nevertheless we conclude that the inclusive distributions are very similar. 
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Figure 20: Comparison of JETSET and our method on an event to event basis: 
We plot here the average of the absolute value of the difference (GeV) between 
the hadronic curves for the same partonic event used in Fig.(|l9|), obtained from 
the two programs as a function of rank, where for every event the particles 
of a certain rank were produced using exactly the same value of z in the two 
alternative procedures. 

But we are definitely not dealing with the same particles on an event to 
event basis. To show this, we compare the two procedures on an exclusive, 
event to event basis in Fig.^. For this plot we devised a procedure to use 
exactly the same values of z for particle production in both programs. Since 
both our procedure, and PYTHIA might need to reject some sequences from 
time to time, we took only those events where no such rejection was required by 
either program. We have plotted the distribution of the absolute value of the 

difference of the two hadronic curves, ~Y^j=\{pc?>) j ~ J2j=i(PPYTHiA)j\ 2 , as 
a function of rank. The curve clearly shows that the hadronic curves obtained 
from the two programs are different. The absolute value prevents cancellation of 
the differences which would occur in an average over a large number of events, 
as we find from the comparison of the inclusive distributions. 



7 Concluding Remarks 

We have presented a precise method to implement the Lund Area Law (Eq. (Eh) 
to fragment a multigluon string state. The final state hadrons are produced in 
an iterative stochastic process. The energy- momentum vectors {pj} can be 
laid out in rank order as a curve, the A-curve, with a vertex vector at every 
point connecting to the directrix curve. The directrix corresponds to the parton 
energy momentum vectors laid out in colour order. It describes on the one hand 
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the orbit of a (massless) g-particle connected to a string and on the other hand 
the whole string surface. 

It is possible to describe the whole process by analytical means, in particular 
in terms of a transfer matrix formalism similar to the one used for the (1 + 1)- 
dimensional model in ||. This time it is necessary to make use of a (1 + 3)- 
dimensional framework. We will investigate this problem in the future. 

In this paper, we have neglected the possibility of introducing transverse 
fluctuations in the hadronisation process. Such fluctuations are introduced on 
the basis of tunnelling arguments in the (1 + l)-dimcnsional model. We feel that 
they are necessary for a consistent quantum mechanical treatment (there is no 
way to localise the string surface area better than what is allowed by Heisen- 
berg's indeterminacy relations). We will examine the effects of such transverse 
fluctuations on our results in the future. 

The directrix curve is in general given by a perturbative cascade (although 
in some cases precise matrix elements are also available). This means that 
the partonic states we are fragmenting, are resolved only down to some cutoff, 
usually in virtuality or in terms of the smallest allowed partonic excitations. 
String dynamics is infrared stable in the sense that minor modifications of the 
directrix only have a local influence, due to the minimal nature of the string 
world surface. We have made use of this freedom in the fragmentation process 
and in this way the partonic states will be defined down to the hadronic mass- 
scale. We will investigate the properties of the soft gluons thus introduced into 
the partonic state, in the future. 

We have shown that the fragmentation process, with step size equal to the 
hadronic mass can be defined also in the limit of a vanishing mass in terms of 
a differential process. The corresponding solution, the V-curve, is stretched in 
between two curves, the original directrix, A, and another curve, the £-curve, 
also with a light-like tangent. The correspondence to the vertex vectors for the 
fragmentation process are connectors to the P-curve reaching out to the A and 
the £-curves. The area in between the P-curve and the directrix A corresponds 
to a generalised rapidity variable in the same way as the average hyperbola 
defines rapidity for the (1 + l)-dimensional model. There is an interesting rela- 
tionship to the group of rotations in a (1 + 4)-dimensional space which we will 
investigate further in the future. 

We have also shown that there is a duality (with properties similar to the 
Parton-Hadron Duality introduced by the St. Petersburg group |l^|) between 
hadronic X-curvc, defined by our process and the original directrix. We note, 
however, that the hadrons produced in our process always stem from the energy- 
momenta of two or more partons. The vertex vectors Xj contain the "memory" 
of the earlier partons. It is nevertheless possible to reconstruct the directrix 
from the X-curve, although the relationship contains a large number of degrees 
of freedom (the number of degrees of freedom increases further if we introduce 
transverse momentum fluctuations in the fragmentation process). We will in- 
vestigate these properties in the future. 

The process we have defined has also been implemented into a Monte Carlo 
simulation program. Due to the way it is constructed it is fairly direct to 
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include also the multi-particle production features of the original JETSET, i.e. 
to include probabilities for different (qq)-fia,vovs, different kinds of mesons, the 
decay of resonances etc. This is necessary for a comparison with experimental 
data. 
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